Kovacs Effect in a Fragile Glass Model 
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The Kovacs protocol, based on the temperature shift experiment originally conceived by A.J. 
Kovacs for glassy polymers Q, is implemented in an exactly solvable dynamical model. This model 
is characterized by interacting fast and slow modes represented respectively by spherical spins and 
harmonic oscillator variables. Due to this fundamental property, the model reproduces the charac- 
teristic non-monotonic evolution known as the "Kovacs effect", observed in polymers, in granular 
materials and models of molecular liquids, when similar experimental protocols are implemented. 
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I. INTRODUCTION 



Systems with glassy dynamics typically exhibit non- 
trivial behavior when they undergo temperature shifts 
within the glassy phase. These systems, being in an out- 
of-equilibrium condition, have properties which are ex- 
pected to depend on their history. This is the 'memory' 
of glassy systems. One memory effect that shows up in a 
one-time observable is the so called "Kovacs effect" [J, 
which manifests itself under a specific experimental pro- 
tocol. This effect has been the subject of a variety of re- 
cent investigations [1, H, ^ IE IS - The characteristic non- 
monotonic evolution of the observable under examination 
(the volume in the original Kovacs' experiment), with 
the other thermodynamic variables held constant, shows 
clearly that a non-equilibrium state of the system can- 
not be fully characterized only by the (time-dependent) 
values of thermodynamic variables, but that further in- 
ner variables are needed to give a full description of the 
non-equilibrium state of the system. The memory in this 
case consists in these internal variables keeping track of 
the history of the system. 

The purpose of this paper is to use a specific model 
for fragile glass to implement the protocol in order to get 
some insight into the Kovacs effect. We show that in spite 
of its simplicity, this model captures the phenomenol- 
ogy of the Kovacs effect, it makes possible to implement 
the Kovacs protocol not only with temperature shifts but 
with magnetic field shifts as well, and allows in specific 
regimes to obtain analytical expressions for the evolution 
of the variable of interest. Furthermore the possibility of 
affording a thermodynamical-like picture through the in- 
troduction of effective parameters can be investigated. 

This paper is organized as follows: in Section II we re- 
view the experimental protocol generating the effect, in 
Sections III and IV we introduce our model and use it 
to implement the protocol, in section V we draw out of 
this model some analytical results, in Section VI an in- 
terpretation of the effect in terms of effective parameters 
is illustrated with final conclusions. An appendix collects 
all terms and coefficients employed in the main text. 



II. KOVACS PROTOCOL 

The experimental protocol, as originally designed by 
A. J. Kovacs in the '60s Jj, consists of three main steps: 

I''* step (1) The system is equilibrated at a given high 
temperature Ti. 

2'"^ step (2) At time t — Q the system is quenched to a 
lower temperature T/, close to or below the glass 
transition temperature, and it is let to evolve a 
period ta- One then follows the evolution of the 
the proper thermodynamic variable (in the original 
Kovacs experiment this was the volume V{t) of a 
sample of polyvinyl acetate, in our model it will be 
the "magnetization" mi[t)). 

3'"'* step (3) After the time ta, the volume, or other cor- 
responding observable, has reached a value equal, 
by definition of ta, to the equilibrium value corre- 
sponding to an intermediate temperature Tf (Ti < 
Tf < Ti), i.e. such that Vr.ita) = ^t this 

time, the bath temperature is switched to Tf. 
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FIG. 1: Kovacs protocol. Starting from an equilibrium con- 
dition ed T — Ti (step (1)) at time t = 0, the system is 
quenched to T = T; and let to evolve (step (2)). In step (3) 
the temperature is switched to Tf . This is done at the time ta 
for which: Vt; (ta) = Vy'. In the frame, the typical evolution 
of the volume V{t) at T = T/, after the temperature switch, 
is illustrated. 
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The pressure (or corresponding variable) is kept 
constant throughout the whole experiment. 

Naively one would expect the observable under consid- 
eration, after the third step, to remain constant, since it 
already has (at time t = f^) its equilibrium value. But 
the system has not equilibrated yet and so the observable 
goes through a non monotonic evolution before relaxing 
back to its equilibrium value, showing a characteristic 
hump whose maximum increases with the magnitude of 
the final jump of temperature Tf — Ti and occurs at a 
time which decreases with increasing Tf — T; . 

We want to implement this protocol on a model for 
both strong and fragile glass first introduced in 0: 
the Harmonic Oscillators-Spherical Spins model (HOSS). 
This model is based on interacting fast and slow modes, 
this property turns out to be necessary for the memory 
effect, object of this paper, to occur. 



III. THE HARMONIC 
OSCILLATOR-SPHERICAL SPIN MODEL 

The HOSS model contains a set of N spins Si locally 
coupled to a set of N harmonic oscillator Xi according to 
the following hamiltonian: 



i=l 



Hxi JxiSi LSi) 



(1) 



The spins have no fixed length but satisfy the spherical 
constraint: X]t=i ~ ^ ■ '^^^ variables are as- 
sumed to relax on a much shorter time scale than the 
harmonic oscillator variables, so the oscillator variables 
are the slow modes and on their dynamical evolution the 
fast spin modes act just as noise. One can then integrate 
out the spin variables to obtain the following effective 
Hamiltonian for the oscillators (for details see 171 , explicit 
expressions of undefined terms appearing in all equations 
hereafter are reported in the Appendix): 
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which depends on the temperature and on the first and 
second moment of the oscillator variables, namely: 
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(3) 



These variables encode the dynamics of the system which 
is implemented through a Monte Carlo parallel update of 
the oscillator variables: 



(4) 



+ n/VN 



The variables are normally distributed with zero mean 
value and variance . The update is accepted according 
to the Metropolis acceptance rule applied to the varia- 
tion Ae of the energy of the oscillator variables, which is 
determined by Ti-eff and, in the limit of large iV, is given 
by: 



Ae KT{mi,7v,2] 



N 



A7712 — Ht (mi, 77i2)A7ni. (5) 



This simple model turns out to have a slow dynamics and 
can be solved analytically. 

Following Q one can derive the dynamical equations 
for mi and 7712 
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(6) 



[/T(mi,m,2) + HT{mi,7n2)7hi\ 



The stationary solutions of these equations coincide with 
the saddle point of the partition function of the whole 
system at equilibrium at temperature T and are given 
by: 



mi 



m2 
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T 
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(7) 



with barred variables from now on indicating their equi- 
librium values. 



Strong and Fragile Glasses with the HOSS 
model 



In spite of its simplicity, the HOSS model allows to de- 
scribe both strong and fragile glasses, characterized re- 
spectively by an Arrhenius or a Vogel-Fulcher law in the 
relaxation time. The following constraint on the config- 
urations space is applied: 



/i2 = m2 ~ mi — mo > 



(8) 



When mg = there exists a single global minimum in the 
configurations space of the oscillators, therefore the role 
of the constraint with mo > is to avoid the existence of 
a "crystalline state" and to introduce a finite transition 
temperature. The stationary solutions for the dynamics 
with this constraint are given by: 



mi 



m2 — mi 



HT{mi,m2) Ht 
Kt {mi, 7712) Kt 

T T ^ 7^ 

KT(fhi,m2) ~ Kt ^ ^ 



(9) 
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T<Tk 



The temperature Tk is determined by: 

Tk = mo Kt^ (771^*= , m^*" ) = mo Kt^ ■ (10) 
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This is the highest temperature at which the constraint 
is fulfiUed, for smaUer temperatures the system relaxes 
to equihbrium configurations which fulfill the constraint. 
For T > Tk therefore the dynamics is not affected by the 
constraint. For T < Tk the system eventually reaches a 
configuration which fulfills the constraint, when this hap- 
pens it gets trapped for ever in such a configuration. This 
is equivalent to have a "Kauzmann-like" transition, oc- 
curring at r = Tfc with vanishing configuration entropy, 
meaning the system gets stuck forever in one single con- 
figuration fulfilling the constraint (see also: Q). 

When there is no constraint, i.e. when mo = 0, then 
Tfc = 0, if the Monte Carlo updates are done with Gaus- 
sian variables with constant variance fi^, this model is 
characterized by an Arrhenius relaxation law: 

Teq ~ (11) 

in so resembling the relaxation properties of strong 
glasses. 

The HOSS model with constraint strictly positive 
(too > 0) can easily be extended to describe fragile 
glasses by further introducing in the variance of the 
Monte Carlo update, the following dependence on the 
dynamics: 



The evolution of the fractional " magnetization" : 



(to2 — to^)(to2 — ml — Too) 



(12) 



In this case the relaxation time turns out to follow the 
generalized Vogel-Fulcher law: 



^Ai/{T-ny 



(13) 



The parameter 7 is introduced to make the best Vogel- 
Fulcher type fit for the relaxation time in experiments, 
making this model valid for a wide range of fragile glasses. 
When the temperature approaches the value Tfc defined 
by I|1U|) . from above, the system relaxes towards config- 
urations close to the ones fulfilling the constraint. The 
variance then tends to diverge, the updates become 
large and so unfavorable, meaning that almost every up- 
date of the oscillator variables is refused. This produces 
the diverging relaxation time following the Vogel-Fulcher 
lawofEq. lO. 



IV. KOVACS EFFECT IN THE HOSS MODEL 

We implement the Kovacs protocol in the model above 
introduced for a fragile glass. The system is prepared at a 
temperature T^ and quenched to a region of temperature 
close to the Tfc, i.e. T/ > Tfc. Solving numerically Eqs. 
© we determine the evolution of the system in both step 
2 and 3 of the protocol. In step 2 the time ta, at which 



(ta) 



is calculated so that: 



- Tf 



(14) 



6mi{t) 



mi{t) 



m^ 



(15) 



after step 3 {t > to) for different values of T/ is reported 
in Figs. 2 and 3 respectively for 7 = 1 and 7 = 2. In all 
the figures, the values for the parameters of the model 
are: J = K— 1, L — H — 0.1, too — 5. For such 
parameter values, the Kauzmann temperature turns out 
to be Tfc = 4.00248. 
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FIG. 2; Fragile glass with 7 = 1. The Kovacs protocol is 
implemented with a quench from temperature Ti — 10 to Ti, 
and final jump (at t — tt) to the intermediate temperature 
Tf = 4.3. The curves, starting from the lowest, refer to Ti = 
4.005,4.05,4.15, the dashed line refers to condition Ti = T/ 
(simple aging with no final temperature shift). 
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FIG. 3; Fragile glass with 7 = 2. The Kovacs protocol im- 
plemented with a quench from temperature Ti = 10 to Ti, 
and final jump (at t = t^) to the intermediate tempera- 
ture Tf — 4.3. The curves, starting from the lowest, refer 
to Ti = 4.005, 4.05, 4.15, 4.25, the dashed line refers to condi- 
tion Ti = Tf (simple aging with no final temperature shift). 



Since the equilibrium value of toi decreases with in- 
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creasing temperature (as opposed to what happens for 
instance with the volume) we observe a reversed 'Kovacs 
hump'. The curves keep the same properties typical of 
the Kovacs effect, the minima occur at a time which de- 
creases and have a depth that increases with increasing 
magnitude of the final switch of temperature. As ex- 
pected, since increasing 7 corresponds to further slowing 
the dynamics, the effect shows on a longer time scale in 
the case 7 = 2 as compared to 7 = 1. 

Actually, since in the last step of the protocol: mi{t — 

to) = and /T/(?Tii,m2) is always positive, from the 
first of Eqs. © one soon realizes that the hump for this 
model can be either positive or negative, depending on 
the sign of the term: 



(16) 



to/ ,m2 



at t = f^. This term is zero when 

so one would expect m2(t = t'^) = rfi^^ to be the bor- 
der value determining the positivity or negativity of the 
hump. Since Ht^ (to^ , 7712) decreases with increasing TO2 

while KTf{fri^^ ,'m,2) increases, it follows that the condi- 
tion for a positive hump is: 



'm-2it = t+) < TO2 



(17) 



For shifts of temperature in a wide range close to the 
transition temperature Tk, where the dynamics is slower 
and the effect is expected to show up significantly on a 
long time scale, the condition m2{t = to) > frt^^ is always 
fulfilled and therefore a negative hump is expected. 



A. Kovacs protocol at constant temperature with 
magnetic field shifts 

Interchanging the roles of T and H , the Kovacs pro- 
tocol can be implemented at constant temperature, by 
changing the magnetic field H instead. From Eqs. @ 
and IjlUI) one can see that the value of the transition tem- 
perature Tk depends on H as well. Different values of H 
determine different values of Tk- Therefore the protocol 
must be implemented in the following way. The temper- 
ature is kept fixed at T^. At this temperature there is 
a "Kauzmann" transition for a specific value of the field 
H = Hk- The temperature Tk decreases with decreasing 
H. So if we work at T = Ti with magnetic fields H < Hk 
we are sure to implement every step of the protocol keep- 
ing the system always above the "Kauzmann" transition 
corresponding to the value of H applied. We start with 
system equilibrated at T = Ti and H = Hi <^ Hk, and at 
time t = we shift instantaneously the field to a larger 
value Hi, such that Hi < Hi < Hk- Then we let the 
system age for a time ta such that TO^'(ia) = ^ . At 
this time the field is shifted to Hf (with Hf < Hi < Hk)- 
The subsequent evolution of the fraction magnetization 



Smiit) is shown in Fig. 4. Again the curves show all the 
typical properties of the Kovacs hump, with a very slow 
relaxation back to equilibrium due to the fact that Hf 
has been chosen very close to Hk- 
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FIG. 4: Fragile glass with 7 = 1. The Kovacs protocol is im- 
plemented at constant temperature Ti = 4.2 with a sequence 
of magnetic field shifts. Hk ~ 2.24787 is the value of the field 
such that Tk =Ti = 4.2. Starting from Hi = 0.1 (step (1)), 
the field is switched to Hi and the system is let to evolve a 
time ta (step (2)). At t = ta, i.e. when mf'{ta) = ^ , the 
field is switched to Hf = 2.17 (step (3)). The curves, starting 
from the lowest, refer to Hi = 2.22, 2.20, 2.18, the dashed 
line refers to condition Hi — Hf- (The values of the other 
parameters of the model are J = K = 1 and L = 0.1) 



V. ANALYTICAL SOLUTION IN THE 
LONG-TIME REGIME 

In the previous Section we have shown, through a 
numerical solution of the dynamics, that the HOSS 
model reproduces the phenomenology of the Kovacs ef- 
fect, showing the same qualitative properties of the Ko- 
vacs hump as obtained in experiments (see for ex. [lllTH') 
or with other models 0,11111111. 

In this section we show that, by carefully choosing the 
working conditions in which the protocol is implemented, 
our model provides with an analytical solution for the 
evolution of the variable of interest. 



A. Auxiliary variables 

In order to ease calculations, as done in H 
convenient to introduce the following variables: 



HT{mi,m2) 

KT{mi,m2) 
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for which the dynamical equations read: 



^2 



- JQt , W2 )/t ("ll , "12 ) 

(1 + DQTimi,m2))fJ.ifTimi,m2) 

2/T(mi,TO2) . 2r f N 

— — + 2^l^fT{ml,m2) 

KT(mi,m2) 



(19) 



We will choose to implement steps 2 and 3 of the proto- 
col in a range of temperature very close to the Kauzmann 
temperature T^. As exhaustively shown in 0,^3 in the 
long time regime the variable fi2it) decays logarithmi- 
cally to its equilibrium value which is small for T ^ Tk- 
So, if ta is very large, the value of the variable ^J.2{t)^ 
which is continuous at the jump, will be small enough to 
fulfill the condition for which the following equation is 
shown to be valid 0: 



d{S^j,2) 



{l + QT{mi,m2)D) 

JQTimi,m2)T 
2(mo -I- fl2) 



(/i2 -I- 6^i2y 

SH2 



(20) 



where now the variable 5^,2{t) — fJ.2{t) — fl2 is used and 
barred variables always refer to equilibrium condition. 
Of course choosing Ti close to Tk and waiting a long 
time ta so that the system approaches equilibrium, allows 
only small temperature shifts for the final step of the 
protocol, meaning that also Tf will be close to Tk . All the 
coefficients which appear in equation (|2()|l (see Appendix 
for complete expressions) in the regime chosen, can be 
assumed constant and equal to their equilibrium values 
with a very good approximation. The equation can then 
be easily integrated to give: 



Hi{5fi2) = exp 



-Ac 



-Cq 



dz exp 



Ac 



2fi(7,7,7+l,-fe) 
2^^1(7,7,7+1, 



7zT 



(21) 



where the superscript + indicates 1 = 1^ and 2F1 the 
hypergeomctric function. This expression simplifies in 
cases 7 = 1, 3/2, 2. All these solutions and relative 
coefficients are reported in the appendix, here we limit 
ourselves to the case 7 = 1 which corresponds to ordinary 
Vogel-Fulcher relaxation law. In this case the solution is: 
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5^12 (t) 
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FIG. 5: Numerical solution (continuous lines) of the Kovacs' 
curves compared to approximate analytical solution at short- 
intermediate (dot-dashed line) and intermediate-long time 
(dashed line). The protocol is applied between Ti = 10 and 
Tj — 4.018. The curves starting from the lowest refer to 
T = 4.005, 4.008. 



One can then expand the variable of interest mi{t) in 
terms of /ii and Sfi2 and obtain the following expression 
for the Kovacs curves: 



6mi(t) = A^ (rn^f ,ml^){iii{t) - fi^) 



(23) 



yl^^(m/', 7712 ('5^2 (<) ~ ^f^t) 



where the coefficients are approximately constant in the 
regime chosen and can be evaluated at equilibrium. 



B. Short and intermediate t — ta 

For small t — ta, a linear approximation for the variable 
Sfi2, with slope given by the second equation of the set 
((T^ evaluated at i = t+, turns out to be very good. 
Inserting this expression in Eq. (|22|l to get /ii(t) and 
then in Eq. H23|) a good approximation of the first part of 
the hump for small and intermediate t — ta is obtained, 
as shown in Fig. 5. 



C. Intermediate and long t — ta 

When t — ta is very large, we can use Eq. I|22l) and the 
pre-asymptotic approximation for /i2(i) (see: 7]) 



fi2it) = (log(iAo) + ^log(log(tAo)) 



-1/7 



(24) 



Inserting this expression in Eq. (|22|l to get /ii(i) and 
then in Eq. H23|l , a good approximation for the hump and 
the tail of the Kovacs curves is obtained. In Fig. 5 we 
show the agreement between the analytical expression so 
obtained and the numerical solution. 



6 




FIG. 6: Effective field vs. effective temperature in the Kovacs' 
protocol. The continuous line AB refers to the last part of step 
2 of the protocol, i.e. aging at Tj = 4.005 after a quench from 
Ti = 10 (we did not present the full line which starts at Te — 
9.13 and He — 0.0826, outside our picture). The continuous 
line CD represents the evolution of the system in step (3) 
of the protocol, after an instantaneous switch of the bath 
temperature from Ti — 4.005 to T/ — 4.018 (resulting in a 
jump from point B to C). The non-monotonicity of the curve 
CD, i.e. of the evolution of He after the jump, is the signature 
of the Kovacs' effect. The dashed line represents simple aging 
at Tf after a quench from Ti. {H = 0.1) 



VI. EFFECTIVE TEMPERATURE AND 
EFFECTIVE FIELD 

The out-of-equlibrium state of the system can be ex- 
pressed by a number of effective parameters which is in 
general equal to the number of independent observables 
considered. In the HOSS model, given the solution of the 
dynamics, a quasi-static approach can be followed to gen- 
eralize the equilibrium thermodynamics (see: 0,^3) by 
computing the partition function of all the macroscopic 
equivalent states at a given time t. The measure on which 
this out of equilibrium partition function is evaluated is 
not the Gibbs measure. One can assume an effective tem- 
perature Te and an effective field He , and substitute the 
equilibrium measure by exTp{—Ti.eff{{xi}, T, He)/Te), and 
in this way determine the partition function and then the 
free energy. The values of the effective parameters at a 
given time t are those that minimize the free energy so 
calculated. In this way one then obtains: 

Te{t) = KT{mi{t),m2{t))[m2{t)-m\{t)] (25) 
He{t) = H - KT{mi{t),m2{t))fiiit) 

We plot He as a function of Te in a Kovacs' setup, in Fig. 
6. We see that in step 2 of the protocol (lower continuous 
line), equivalent to a simple aging experiment, the effec- 
tive magnetic field relaxes monotonically to the value H. 
In step 3 (upper continuous line CD), after the final jump 
of the bath temperature, represented in the figure by the 
jump from point B to point C , the effective magnetic 
field goes through a non-monotonic evolution before re- 



laxing to the equilibrium value H. This is when the "Ko- 
vacs" effect occurs. A conclusion that can be drawn is 
that a thermodynamic-like picture in terms of only the ef- 
fective temperature is not possible in the Kovacs setup if 
not at cost of neglecting effects of the order of magnitude 
of the "Kovacs effect" itself. So, while in an aging exper- 
iment in the long time regime He(t) — H is very small 
compared to Te (t) — T (so that one can consider He = H 
and use only Te as effective parameter), in the Kovacs 
protocol, it is in the non-monotonic evolution of the ef- 
fective field that the memory effect manifests itself. An 
additional effective field is then needed to recover a com- 
plete thermodynamic-like picture of the system inclusive 
of the Kovacs' effect. The dashed line in the figure rep- 
resents a simple aging experiment at T = Ty, in this case 
a thermodynamic-like picture with only an effective tem- 
perature is possible, assuming He ~ H. One can argue 
from the figure that such a picture would be possible also 
in step 3 of the protocol (curve CD) since the two curves, 
for Te close enough to Tf, coincide. But this happens 
when the system is close to equilibrium and the signa- 
ture of the memory effect, the non-monotonic evolution, 
is lost. These conclusions confirm the results obtained in 
Ref. 13 , where the impossibility of a thermodynamic-like 
picture with only the effective temperature was based on 
a potential energy landscape (PEL) analysis. The molec- 
ular liquid studied in |6j , in the last step of the Kovacs 
protocol, explores regions of the PEL never explored in 
equilibrium, and so a simple mapping to an equilibrium 
condition at a different temperature (the effective tem- 
perature by definition) is not possible. 



VII. CONCLUSION 

We have shown that a simple mode with constrained 
dynamics like the HOSS model, is rich enough to repro- 
duce the Kovacs memory effect, even allowing to obtain 
analytical expression for the Kovacs hump in a long time 
regime. The Kovacs effect is observed in many experi- 
ments and models, showing common qualitative proper- 
ties which we have found to be shared also by the model 
analyzed in this paper. The quantitative properties de- 
pend on the particular system or model analyzed. 

As far as it concerns the HOSS model, it turns out 
that for the slow modes, i.e. the oscillator variables, 
fixing the overall average value, the magnetization mi, 
does not prevent the existence of memory encoded in 
the variable 7712 , which keeps track of the history of the 
system. The equilibrium value of m2 increases with tem- 
perature while the equilibrium value of toi decreases with 
increasing temperature. Therefore after the final switch 
of temperature, since m2(ta) > w^^, the variable TO2 has 
a value corresponding to an equilibrium condition at a 
higher temperature (memory of the initial state at tem- 
perature Ti) so driving the system towards a condition 
corresponding to a higher temperature, i.e. smaller val- 
ues of mi, determining the hump. 



It is important to stress that a fundamental ingredient 
in the HOSS model is the interaction between slow and 
fast modes. Due to this interaction the equilibrium 
configurations of the oscillator variables at a given 
temperature are determined by both m2 and mi, the 
first and second moment of their distribution, whose 
dynamical evolution is interdependent. When such 
interaction is turned off (by setting J = ) essentially 
only one variable is sufficient to describe the equilibrium 
configurations and the dynamics of the system, and 
the memory effect is lost. In this respect this model 
constitutes an improvement to the so-called oscillator 
model |l2j| within which such memory effect cannot be 
reproduced. Another important conclusion, confirming 
previous results P| , which can be drawn from this model 
is that a complete thermodyanmic-like picture inclusive 
of the Kovacs' effect, with only an effective temperature 
is not possible and that also an effective field in this 
case is needed. In the present model one can also study 
temperature cycle experiments of the type carried out in 
spin glasses (see. [T^'l. leaving room for further research. 

G.A. and L.L. gratefully acknowledge the European 
network DYGLAGEMEM for financial support. 
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^ T{c + n) n\ 



Bl = exp[AQ 



2^1(7,7,7+1,-^)^ 



7('5/i2^)^ 

J Qtj {ml' , fhl' )Tf J Qtj Tf 



2{mo + fi2) 2(mo + M2) 

{WT + T/2)Kt 
mi(Jmi +L+{wT + T/2)Kt) 
2miAip{mi, 7712) 
I + DQt 
87 1 + DQt + Pt 



In this Appendix we report all the explicit expressions 
for terms appearing in the text. In Eqs. (O and © we 
have: 



1. Solutions of Eq.lHoll for 7 = | 



Wt (toi , ^2 ) ^ y/J^m2 + 2JLmi + + T^/i 

J2 



KT{mi,m2) = K 



HT{mi,m2) = H 



WT{mi,m2) +T/2 
JL 

WT{mi,m2) + T/2 



/T(mi,TO2) 



a^KT{mi,m2) 
2f 



Erfc [aT{mi,m2)] 



/T(mi,TO2) 



■ exp [ay(TOi, m2) — a|n(mi, TO2)] 

(T^KT{mi,m2) ^ r ^ I M 
Erfc [aT{mi,m2)\ 



KT{mi,m2)wT{mi,m2) /t(»71i,TO2) 



where: 



~ t N 2 , ,HT{mi,m2) 2 

WT(mi,m2) = m2-mi + {—— ™ 

KT{mi,m2) 



aT{mi,m2) = J — — r 

y 8wT[mi,m2) 

aT{mi,m2) _ 2KT{mi,m2)wT{mi,m2) 
aT{mi,m2) T 

In Eqs. mi), C21), 1201, (ED, (123 and 



- 1 



Aii('5Ai2) = 



1- a/1 



61.12 



1+x/l 



2Aq 



S^^2 



1 + Jl 



_3/2 
*'2 



-Co / dz 

Is^i \I-./I 



\ 



I 



2. Solutions of Eq.(|20ll for 7 = 2 



f^l{SH2) 



SfJ.2 



5^J.2 + fJ.2 



-Co I dz 











2 \ "2 


e '^i'l+A 
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